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An introductory exposition of Chandrasekhar's gravitational dynamical friction, appropriate for 
an undergraduate class in mechanics, is presented. This friction results when a massive particle 
moving through a "sea" of much lighter star particles experiences a retarding force du to an ex- 
change of energy and momentum. General features of dynamical friction are presented, both in an 
elementary and in a more elaborate way using hyperbolic two-body interactions. The orbital decay 
of a massive particle in an homogeneous gravitational system is solved analytically, that leads to 
an underdamped harmonic oscillator type of motion. A numerical integration of the equation of 
motion in a more realistic case is done. These results are compared to those of an iV-body computer 
simulation. Several problems and projects are suggested to students for further study. 



I. INTRODUCTION 

Classical mechanics, perhaps the oldest of the phys- 
ical sciences, continues to be an area of intensive re- 
search, both in its foundations^* 2 " and applications^ 3 ^ and 
a source of discussion and examples in teaching. Ap- 
plications range from the modeling of cellular mechan- 
ical processes^ to solar system dynamics 3 and galactic 
systems £L£ 

In describing nature students learn from their first 
courses, and particularly in laboratory experiments, that 
"the forces on a single thing already involve approxima- 
tion, and if we have a system of discourse about the real 
world, then that system, at least for the present day, must 
involve approximations of some kind" ; as mentioned by 
Feynman on introducing the subject of friction^ 

This phenomenon is usually introduced in text- 
books 10,11,12,13 and lectures by considering the slide of 
a material block on a surface, and a distinction between 
static and kinetic friction is made. A classical example of 
the effect of a friction-like force is the motion of a mass 
attached to a spring inside a viscous medium, where the 
corresponding differential equation is solved, and its be- 
havior studied. At the end, one invariably needs to state 
that friction and its origin is a complicated matter, in- 
volving complex interactions at the atomic and molecular 
level among the surfaces in contact (i^iSiiSiii 

Several non-typical examples of mechanical friction for 
introductory courses exist f i2*i2i 2 iii2I that help both teach- 
ers and students alike in lectures on mechanics. All fric- 
tion related problems are a background for discussing the 
important connection between the work-energy theorem 
and dissipative systemsi 22 i 23 i 24 

The purpose of this paper is to bring an example from 
astronomjiSSiSLSl closely related to standard mechani- 
cal friction, namely: dynamical friction. This process 
was first introduced in stellar systems by Subrahmanyan 
ChandrasekhariS&Si In brief, a massive particle m experi- 
ences a drag force when moving in a "sea" of much lighter 
star particles m* by exchanging energy and momentum. 



An elementary understanding requires only some basic 
ideas from mechanics, and hence suitable for presenta- 
tion in introductory courses. 

Dynamical friction is important in astronomical stud- 
ies of, for example: the fate of galaxy satellites 30,31,32 or 
globular clusters^ 3 - orbiting their host galaxies, the sub- 
structure of dark halos surrounding galaxiesj 3 ^ 3 ^ 3 ^ and 
the motion of black holes in the centers of galaxies. 3 ? 
It has been proposed to explain the formation of bina- 
ries in the Kuiper-behjiS, and the migration of Jupiter- 
mass planets in other solar systems from the outer parts 
where they presumably formed (> 1 AU) to the small or- 
bital distances (< 0.1 AU) at which they are observed* 3 ^ 
It even has been considered in the motion of cosmic 
strings 4°» 

The presentation of this topic to students, in a lower 
or upper-undergraduate class on mechanics^i*^ or com- 
putational physics^ 3, will enhance their appreciation of 
physics in describing nature and expose them to another 
example of classical mechanics. Furthermore, students 
will obtain a glimpse of an area of astronomical research 
important for the understanding of the fate and behavior 
of stellar systems. 

The organization of this paper is as follows. In Sec- 
tion II basic elements of the theory of dynamical friction 
are presented. Firstly, elementary arguments are used 
to elucidate them. Secondly, Chandrasekhar's approxi- 
mation using two-body hyperbolic Keplerian collisions is 
considered. In Section III a simple analytical problem for 
the motion of a massive particle in an ideal homogeneous 
stellar system is solved; a damped harmonic oscillator is 
found. In Section IV a more realistic astronomical exam- 
ple that requires the numerical integration of the equa- 
tion of motion is presented. Comparison with a computer 
experiment is done afterwards. Final comments as well 
as some ideas for problems and projects of further study 
are provided in Section IV. An appendix contains some 
astronomical units and standard units used in gravita- 
tional computer simulations. 
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FIG. 1: Illustration of the deceleration a heavy particle m ex- 
periences when moving in an homogeneous and infinite "sea" 
of much lighter particles to*. 



II. DYNAMICAL FRICTION 

Two equivalent approaches to compute the dynami- 
cal friction a massive particle to experiences as it moves 
through a stellar system of much lighter stars to* are the 
following. 6 (1) Particle to produces a region of star over- 
density behind it, much like the wake behind the motion 
of a ship, that in turn exerts a gravitational pull on m 
leading to its deceleration. 44 (2) Particle to moves in the 
"sea" of lighter particles to* and an energy exchange oc- 
curs, increasing that of the lighter ones at the expense of 
the heavy one leading to a breaking force for to. In the 
latter picture the basic features of dynamical friction are 
easier to compute and understand by elementary meth- 
ods. Here the latter picture is taken. 



A. Elementary estimate 

Consider a particle to moving with velocity v m in an 
homogeneous background of stationary lighter particles 
of equal mass to* ; see Figure ^ Assume only changes 
in kinetic energy. As to moves through, a particle to* 
incoming with impact parameter b will be given a velocity 
impulse of about the acceleration a times the duration of 
the encounter At. This can be approximated as 



Av* 



Gm 



The kinetic energy gain of to* is therefore 

1 ,a ^2 1 (Gm\ 
Afi» w -m*(At>*J w -m* I — — I 



(1) 



(2) 



The total change in velocity of the massive particle is 
given by accounting for all the encounters it suffers with 
particles to* . The number of encounters with impact pa- 
rameter between b and b+Ab is AN ss no (v m At) A(nb 2 ): 
where hq is the number density of background stars. The 
total change in velocity of to at the expense of the energy 
lost by stars is then 



1 



dt 



cLE* 



-dN 



irG 2 pom 



db 

T 



(3) 



where we set po = riom*, the background density, and 
frmin an d frmax &te a minimum and maximum impact pa- 
rameter, respectively. Letting In A be the resulting inte- 
gral, the deceleration of to due to its interaction with an 
homogenous background of particles stars is 



dv m irG 2 pQin 
dt vL 



In A. 



(4) 



The velocity impulse on to* has a perpendicular At>j_ 
and parallel Av\\ component; see Figure^ It is not diffi- 
cult to see that a mean vector sum of all the Av± contri- 
butions vanishes in this case. This is not true however for 
the mean square of Av± £^ Thus the dynamical friction 
force is along the line of motion of to. 

Several key features of dynamical friction are observed 
from equation Q in this elementary calculation, that 
appear also in more elaborate treatments. (1) The de- 
celeration of the massive particle is proportional to its 
mass to, so the frictional force it experiences is directly 
proportional to m 2 . (2) The deceleration is inversely pro- 
portional to the square of its velocity v m . 



B. Chandrasekhar formula 

A further step in calculating the effect of dynamical 
friction is to consider hyperbolic Keplerian two-body en- 
counters. Such analysis was done by Chandrasekhar 
The resulting formula is provided in textbooks on stel- 
lar dynamics^ For completeness such calculations is pro- 
vided here, following Binney & Tremaine. 

Use of well known results from the Kepler problem for 
two bodies in hyperbolic encounters are usediiS^ii* 1 ^ 1 ^ 4 ^ 
The two-body problem can be reduced to that of the 
motion of a particle of reduced mass p — toto*/(to + to*) 
about a fixed center of force: 



pr 



(5) 



where k = Gram*, r = r* — r m is the relative vector 
position of particles to and to* , and f its unit vector; see 
Figure |21 The relative velocity is then V = v* — v m , and 
a change in it is 



AV = Av, - A^ 



(6) 



The velocity of the center-of-mass of to and to* does 
not change, hence 

to* Av* + to Av m = . (7) 

From equations JHJ) and {JJ the change in velocity of to 



Av m = -( 



TO + TO* ' 



AV. 



(8) 



Once AV is determined, Av m can be found from equa- 
tion iJSJ. From the symmetry of the problem, is better 
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FIG. 2: Dispersion of a "reduced" mass particle /i in the 
potential of a static body of mass [m + m») moving with 
relative speed Vq. The scattering angle is 9. 



to decompose AV in terms of perpendicular and parallel 
components: 

AV = AV,,+AV ± , (9) 



with 



AVj| | = V a cos 6 and | AV± \ = V sm9, (10) 



where 9 is the angle of dispersion and Vq the initial speed 
at infinity; this being the same after the encounter since 
only kinetic energy changes are considered; see Figure |21 
From geometry, the angle a in Figure is related to the 
orbit's eccentricity e by^& 



1 



cos a 



cot • 



e 2 
where 9 + 2a — it. Physically e is given by 



e = k 1 



2EL 2 



(11) 



(12) 



where E = fiV 2 /2 in the kinetic energy and L = [ibVo 
the angular momentum magnitude. Since 



sin 9 



2 tan | 1 - tan 2 % 

jr-n and cos 9 = —f- 

1 + tan 2 I 1 + tan 2 | 



after some algebra it is found that 
2bV 3 



AV, 



G(m + m») 



1 



G 2 (m + m*) 2 



|AF H | = 2F 



1 



. b 2 V± 
G 2 (m + m*) 2 



(13) 



(14) 



(15) 



Using equation © the perpendicular and parallel mag- 
nitudes of the components of Av m follow: 



|Av r , 



2bmV 3 
G(m + m») 2 



b 2 V* 



2mV 

l AVm IM = 7 TTTT 

(m + m*) 



G 2 (m + m*) 2 



b 2 V± 



G 2 (m + m*) : 



, (16) 



(17) 



In a homogeneous sea of stellar masses all perpendicular 
deflections cancel by symmetry. However, the parallel ve- 
locity changes are added and the mass m will experience 
a deceleration. 

The calculation of the total drag force due to a set 
of particles m* is as follows. Let /(v*) be the number 
density of stars. The rate at which particle m encounters 
stars with impact parameter between b and b + d&, and 
velocities between v* and v» + dv* , is 



27r6d6-F -/(v»)d 3 v, 



(18) 



where d 3 v„ is the volume element in velocity space. The 
total change in velocity of m is found by adding all the 
contributions of | Av m || | due to particles with impact pa- 
rameters from to a 6 m ax and then summing over all 
velocities of stars. At a particular v* the change is 



dv l; 



dt 



= V -/(v„)d d v, 

v. Jo 



|Av m ,||27r6d6. (19) 



The required integral is 



I 



1 



b 2 V£ 



-i -l 



1 



G 2 (m + m*) 2 
bdb 1 
2a 



ib 2 



bdb 

ds 

s 



where a = Vq /G 2 (m + m 1 ,) 2 and s = 1 + ab 2 , with s n 
1 + a6 2 nax . Evaluating the integral yields 



lG 2 (m + m») 2 . r 
1 = -. In 1 

2 V n 



A 



where 



A 



V 2 



G(m + m*) b 

min 

Putting these results together in equation (|19|l 

<lV "' 27rG 2 ln(l + A 2 )m,(m + m,) 

x /(v^d^v*- 13 . 

v* - v m \ 6 



(20) 



dt 



(21) 



The quantity In A is called the Coulomb logarithm in 
analogy to an equivalent logarithm found in the theory of 
plasma. The factor In A reflects the fact that the cumu- 
lative effect of small deflections is more important than 
strong or close encounters. This may be seen geometri- 
cally from Figure |3 were the stronger the deflection the 
smaller is the parallel component contributing to the slow 
down of m. 

The determination of the limits 6 min and 6 max is not 
an easy matter and depends on the problem at hand. In 
this approximation 6 m ; n satisfies Vq = Gm/b m i n , where 
Vq depends on the relative velocity of m and m„. If the 
motion of m is relatively slow in comparison to that of 



the stars, Vq can be approximated for example by the 
root-mean-square value velocity of stars V rms . The outer 
limit & m ax is in principle the radius at which stars no 
longer can exchange momentum with m. If m is close 
to the center of a stellar system & max can be taken as a 
particular scale-radius of the system; for example, where 
the star density falls to half of its central value. 

In typical astronomical applications A> 1. For ex- 
ample, consider the motion of a massive black hole of 
mass m w 10 5 M Q near the center of a dwarf galaxy. 
These galaxies have V rnls w Vq « 30 km s _1 , character- 
istic radii 6 max ~ 3 kpc and stars of masses m* ps 1M©. 
Using these values we obtain A ps 6.3 x 10 3 . This allows 
to use the approximation ln(l + A 2 ) ps 2 In A. Note that 
In A shows a weak dependence on Vq that is usually ne- 
glected. Values of 2 < In A < 20 are typically found in 
astronomical literature. 

Now, the integration of equation (|21|l over the velocity 
space of stars is required. Writing equation l|21l) as 



dv m 
dt 

p(v„; 



= G 



p(v*)(v» - v m ) 



(22) 



4nG(m + m*)m* In A/(v*) 



it is noticed that represents the equivalent problem of 
finding the gravitational field (acceleration) at the "spa- 
tial" point v m generated by the "mass density" p(v*). 
From gravitational potential theory^S^i the acceleration 
at a particular spatial point r is given by 



a(r) = G 



p(r')(r' — r) d 3 i 



Gr 



/ p(r') dV . 
Jo 



This is the known result that only matter inside a par- 
ticular radius contributes to the force. In analogy to the 
gravitational case, the acceleration is given by the total 
"mass" inside i>* < v m , is 



dv„, 
dt 



Gv„ 



p(v*)d 3 v* 



For an isotropic velocity distribution: 

dVm 



dt 



Cdf 



-Cdf v„ 



(23) 



167r 2 G 2 TO*(m + m*) In A / f(v*)v*dv 



This is called Chandrasekhar dynamical friction formula. 
It shows that only stars moving slower than v m con- 
tribute to the drag force on the massive particle. 

If stars have a Maxwellian velocity distribution func- 
tion, 



/(«.) = 



no 



(2trj 2 ) 3 / 2 



(24) 



the integral in i|23|) in done by an elementary method. In 
dimcnsionlcss form it is 




FIG. 3: Forces acting on a massive particle m moving with 
velocity v inside a stellar system: F g is the gravitational force 
and Fdf the dynamical friction force. The latter acting oppo- 
site to the direction of motion of m. 



where y 2 — v 2 /2o 2 and X 
parts results in 



J(y/2a). Integrating by 



T - n ° 
m " 47T 



Erf (X) - H e~ x2 

\/7T 



where Erf(x) = (2/y/n) Jq e~ v dy is the error function. 
If po = nom«, the density of the background of stars, and 
assume that m 3> m*, the deceleration of m inside an 
homogeneous stellar system with isotropic velocity dis- 
tribution is: 



dt 

r df 



-r df v 

47rG 2 In Ap m 



(25) 



Eri(X)-^=e- x2 



III. AN ANALYTICAL EXAMPLE 

A simple application of Chandrasekhar's formula l|25|l 
for an homogeneous spherically symmetric stellar system, 
although not infinite, is presented. The problem consists 
in determining the motion of a massive particle m subject 
to gravitational and dynamical friction forces. The stellar 
system has a radius R and total mass M; see Figure |3] 
The equation of motion for m is 



d 2 !- 

'dP 



'df 



-rriVtp(r) + madt , 



(26) 



where tp(r) is the gravitational potential, and adf is given 
by equation (J5SJ. 

Equation (|26|) is not in general tractable by analytical 
methods, so some approximations are required. Zhac 35 
has found an approximation to the term associated with 
the velocity distribution in equation (|25|l . namely: 



5 



that works to within 10 percent for < X < oo. When 
to moves slow in comparison to the velocity dispersion of 
stars, v <C a, x(0) ~ 3/4, and when v — a/2 x(l) w 3/7. 
Note that in the case of a very fast relative motion of m 
dynamical friction is negligible; a situation analogous to 
when a block of material slides fast. Using the previous 
approximation, and considering x = 3/4, the frictional 
force Fdf in equation i|26|) becomes: 



'df 



3nG 2 \nAp m 2 



(V2a) s 



■ v = -7V, 



(28) 



where 7 = 3ttG 2 In Ap Q m 2 / (y/2a) 3 . 

To determine F g recall that the potential is related to 
the density through Poisson equation, 



V 2 c/?(r) = AnGp{r) 



(29) 



whose solution for a spherically symmetrical system of 
radius R is 



(p(r) = -AirG 



p(r)r 2 dr 



p(r)r dr 



In a constant density po system the potential is 

<p(r) = -2kG P , (R 2 - V) , 

and the gravitational force on to is 

4 



-rriV(p(r) 



-ixGpam r = — k r . 



(30) 



(31) 



(32) 



with fc = AnGpom/S. This is the well known result from 
introductory mechanics that a particle inside an homoge- 
neous gravitational system performs a harmonic motion. 

Combining equations QSt an d <|3"2*f the resulting equa- 
tion of motion is 



d 2 r , 
to-— + kr 

dt z 



7 v = . 



(33) 



This is the same equation, for example, as that of a 
mass attached to a spring with stiffness constant k in- 
side a medium of viscosity 7; i.e., a damped harmonic 
oscillatoriifiiiiii 2 *^ The solution of equation (|33(l in a 
plane, under arbitrary initial conditions 

x(0) = x Q , x(0) = u ; y(0) = y , y(0) = v Q , 

where the dot indicates a time derivative, is^Sd^ 

e -(/3+S)< 



x(t) 



[0 



2R 

2?t . 



l)^+(e 2 «* + l)e] x Q } 



-(0+Ot 



2^ 



{0 



l)v - 



where 2/3 = 7/m and ^ = \J (3 2 — w 2 , with w 2 = fc/m. 

The behavior of m is dictated by the relative values of /? 
and ui. The values of In A and a are first to be estimated. 
Take & max ~ R and & m ; n = Gtu/Vq. An estimate of Vo 
may be obtained from the virial theorem^ that relates 
the kinetic T and potential energy W of the system by: 



2T = -W . 



(35) 



For an homogeneous system of size R the potential energy 
is 



r R 3 GM 2 

W = ~4nG p Mrdr = --—^-, (36) 
Jo 



5 R 



and the kinetic energy is taken as T = MV 2 /2. This 
leads to 



Vn 



3 GM 



5 R 



3a z 



(37) 



where the last term provides an estimate of the one- 
dimensional velocity dispersion under the assumption of 
isotropy in the velocity distribution of stars. Using equa- 
tion (|3T|) f3 and u are: 



P 



45 [5 G 2 mM In A 



16 V 2 R 3 (GM/R) 3 / 2 ' 



w 



GM 



(38) 



The resulting Coulomb logarithm is In A = ln[3M/(5m)]. 

To compare the numerical values of (3 and coq is better 
to use another system of units than a physical one. Let 
G = M = R = 1, that is a common choice in A^-body 
simulations in astronomy; to return to physical units one 
can use Newton's law and set G to the appropriate value 
(see Appendix). In these units, relations (|38(l become 



R 45 5 1 ( 3 

16 V 2 V 5m 



lj = 1. 



(39) 



If m = 1/100 then In A w 4 and /? w 0.2 < w - Hence an 
underdamped harmonic motion for the massive particle 
results. If to = 1/10 then In A ps 2 and /3 » 0.8 » Wq, 
so the motion of to will be strongly damped. Note that 
an upper limit to to is set when to = 3/5, leading to 
In A = 0; i.e., no dynamical friction results. For larger to 
a negative f3 is obtained. Clearly, the model fails and the 
behavior of the dynamics is unrealistic. 

For cases of interest, where m <§C M, it follows that 
f3 < bjQ and the resulting motion (|34|l , after some algebra, 
is 



x(t) 

y(t) 

where to 2 



u + (3x . ' 
xq cos ojl H sin tot 

, v Q + (3y . 
yo cos cot H sin cut 



-Pt 



(40) 



f3 2 . Note that a time-scale when the 



[(e 2?t - l)[3 + (e 2? * + 1)^] y } , (34) orbit decays 1/e is given by r df = 1/(3. 
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(y ,x ) = (0.59,0.44) 
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FIG. 4: Orbital decay of a massive particle m in a homoge- 
neous stellar system due to dynamical friction. Initial con- 
ditions are (xo,yo) — (0,0.59) and (uo,vo) = (0.44,0). Left 
panels are for m — 0.01 and right ones for m — 0.02. Top 
panels show the orbit (solid line) in the xy-plane and bottom 
ones the time evolution of the distance r from the center. The 
dynamical friction time scale Tdf = 1//3 is indicated by an ar- 
row. Dashed lines are orbits without considering dynamical 
friction. 



Take as example m — 0.01. If the initial position of 
m is at (x, y) = (0.59, 0.0) and its velocity (0.0, 0.44) the 
resulting orbit is that shown in Figure QJZe/it with solid 
line. Doubling m results in the orbit shown in Figure 0} 
right. The dashed lines in correspond to the orbit of m 
without dynamical friction. The effect of increasing the 
mass of m on the orbit, and on the decay time Tdf, is 
clearly appreciated. 

This example shows the basic features of, for example, 
the orbital decay of a satellite galaxy toward the cen- 
ter of its host larger galaxy. It may be applied also to 
the motion of a massive black hole near the center of a 
galaxy or star cluster, where to some approximation the 
gravitational potential can be taken as harmonic. More 
realistic situations require however the numerical integra- 
tion of the orbit and/or an A-body computer simulation. 
A particular case of these are treated next. 



IV. A MORE REALISTIC EXAMPLE 

Chandrasekhar's formula H25|) although derived as- 
suming an infinite homogeneous system may be ap- 
plied, to some degree, when stellar systems are non- 
homogeneouswi In this case, local values for the density 
p(r) and the velocity dispersion <r(r) are used. Here the 
motion of a massive particle m inside a non-homogeneous 



stellar system is considered, both using a semi-analytical 
method and A-body simulation, to illustrate further the 
application of dynamical friction. 



A. Semi-analytic treatment 

A simple representation of a stellar system, such as a 
globular cluster or an elliptical galaxy, is provided by the 
Plummcr model. Its potential and stellar density are, 
respectively**^ 



GM 



( r 2 +a 2)l/2> 



p(r) 



3Ma 2 /4ir 

( r 2 +a 2)5/2- 



(41) 



where M is the total mass, and a the scale-radius of 
the system. In a spherical system with isotropic velocity 
distribution the equation of "hydrostatic" equilibrium 62 
is satisfied: 



ld( P( 7 2 ) 
P 



dr 



dip 
dr 



u ^ 2 / \ 

3T -» o- ( r ) 



6 



(42) 



The last result follows from noticing that p oc p 5 , and 
imposing boundary conditions that both pa 2 and p go 
to zero at infinity. 

Equations (|4"T)l and j3H will be used in equation lt231) 
to compute the orbital motion of a massive particle m. It 
rests to determine 6 m i n and 6 max . The former is evaluated 
at local values, & m j n = Gm/[3a 2 (r)], and the latter is set 
fix to 6 max = a- 

The equation of motion (|26|l for m can now be inte- 
grated numerically using standard methods^S^i or using 
the one discussed by Feynman^ for planetary orbits (§9). 
Here a fourth-order Runge-Kutta algorithm with adap- 
tive time-step was used. The initial conditions for m are 
the same as those used in the analytical case. 

In Figure [S] the resulting orbit from the numerical in- 
tegration is shown as a dashed line. Also, the behavior of 
the x and y coordinates, and of the distance r of m to the 
center, as a function of time are shown. The typical de- 
cay of the orbit is evident. In the same figure results from 
an iV-body simulation are displayed, that are described 
next. 



B. N-body simulation 

The use of iV-body simulations allows to study more 
realistically the different dynamical phenomena that oc- 
cur in stellar systemsA^ Several iV-body codes with dif- 
ferent degrees of sophistication have been developed for 
astronomical problems in mindi^S^iS Some low- A'' sim- 
ulations can be run nowadays using a personal computer 
with publicly available A-body codesiS 

Barnes' tree-code in Fortran, and some of his pub- 
lic subroutines are used to simulate the motion of m 
inside a Plummer model. A numerical realization of 
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FIG. 5: iV-body simulation of the orbital decay (solid line) 
of a massive particle m inside a Plummer stellar model. The 
semi-analytical (dashed line) and the analytical calculation 
(dotted line) of Figure 21 are drawn for comparison. These 
overestimate the effect of dynamical friction in comparison to 
the numerical simulation. 



this model with N — 10 5 particles is used with indi- 
vidual "star" masses of m* = X/N. The massive particle 
m = 1/100 with initial conditions (yo,Xo) = (0.59,0.44) 
is set "by hand" inside the numerical Plummer model. 
In iV-body units the scale radius is a — 3ir/lQ = 0.59 
(see Appendix). 

The circular period at radius r is r = 2Trr/V c (r), where 
the circular velocity and integrated mass for a Plummer 
model are given, respectively, by: 



V c (r) 



GM(r) 



M(r) 



M(r/a) 



[1 + (r/a) 2 ] 3 / 2 



From this, an orbital period of r a = 4.8 time units at 
r = a results. The simulation was run for t = 10 ~ 2r a 
time units. The parameters for running the tree-code 
in serial were those provided by Barnes at his Internet 
site for an isolated Plummer evolution. The quadrupole 
moment in the gravitational potential is activated. The 
simulation took about 5.2 CPU hours on a PC with an 
Athlon 2. 2 GHz processor, and 512 KB of cache size. En- 
ergy conservation was < 0.04 percent, that is considered 
very good. 

Figure |S] shows the orbital evolution of the massive 
particle m in the A^-body system as a solid line. The 
dashed line corresponds to the semi-analytical calculation 
of Section llV Al This follows closely the orbit of m in the 
Af-body simulation for about r a time units. Afterwards, 
it deviates from the A^-body result. In the r-t panel, the 
analytical solution l|40|l is shown as a dotted line; that is, 
assuming the total system was homogeneous. 



Both approximations overestimate the decay rate of m 
in comparison to the N-ho&y simulation. Taking 6 max = 
a/5 leads to a somewhat better agreement, but does not 
reproduce the A^-body result. Rather surprisingly, the 
analytical result does a fair job in reproducing the overall 
orbital decay in this case. 



V. FINAL COMMENTS 

The approximations in deriving Chandrasekhar for- 
mula limits, obviously, its application to more complex 
stellar systems than the one considered here. However, 
it is remarkable that equation l|21|l leads to reasonably 
well results when used with values under a local approx- 
imation. 

In similar vain to the study of the friction between 
surfaces ^i&ii dynamical friction is a complex subject. 
Elaborate calculations based on Brownian motion^ lin- 
ear response theory, resonances, and the fluctuation- 
dissipation theorem exist i2S*SS These that are steps for- 
ward toward a more complete physical theory for this 
process. 

Instead of listing explicitly some of the shortcomings 
of Chandrasekhar dynamical friction formula 6 when ap- 
plied to gravitational systems, the student is encourage 
to think on some of them and possible improvements on 
such formula. 

From the point of view of an introductory or interme- 
diate class on mechanics the exposure of students to non- 
typical problems, as the one presented here contributes 
to further their understanding and appreciation of the 
subject. 

Some ideas that may lead to problems and/or projects 
for students are: 

1. How would the analytical solution considered here 
would be changed if the Plummer model is used? 
What type of approximations would be required to 
make? How does In A change? 

2. If <j is a measure of the kinetic energy per unit 
mass of stars, what is an estimate for its mean in- 
crease due to the energy lost by the massive particle 
during its decay? 

3. How would the orbital decay time be changed for 
different types of initial eccentricities of the massive 
particle? 

4. Consider a star cluster (ra = 10 6 M©) in circular 
orbit at a distance of r = 5 kpc from the center of 
our galaxy (M « 6xlO n M Q , R » 150 kpc). Would 
it be expected to fall to the center within the age 
of the universe, say t = 10 10 yr? Typical velocities 
for stars and dark matter particles at that distance 
are about 200 km/s, and the scale-radius may be 
around 5 kpc. What if instead of a star cluster 
we have a galaxy satellite, such as the Magellanic 
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TABLE I: Astronomical units 



Unit 



Astronomical unit" 
Parsec 

Kiloparsec 
Solar mass 
Year 



"Mean sun-earth distance 



Equivalence 



AU = 1.496 x 10 11 m 
pc= 2.063 x 10 5 AU 

= 3.261 light-years 
kpc=10 3 pc 
M Q = 1.989 x 10 30 kg 
yr = 3.156 x 10 7 s 



TABLE II: From A r -body units to astronomical 



Stellar system 


Ul 




u v 


u t 






M 


km/s 


Myr 


Globular cluster 


50 pc 


10 b 


9.3 


5.3 


Galaxy 


10 kpc 


10 11 


207.4 


47.2 


Cluster of galaxies 


5 Mpc 


10 15 


927.4 


5271.4 


The Gravitational 


constant 


can be 


expressed 


in terms 



of typical astronomical values, for example, as: 



Clouds, with m 
100 kpc? 



10 10 Mm and at a distance of 



5. How do results change if instead of a Plummer 
model a more pronounced density profile is used, 
such as the HernquisliS model? How does the num- 
ber of particles TV in a simulation affect the decay 
rate? 

6. As the massive particle moves through the stellar 
system it induces a density wake behind it. Can 
this be detected in an iV-body simulation on a home 
computer? How about looking for this wake in the 
phase-space diagram (e.g. a plot of x-x) of stars 
near the the massive particle? 

7. How good do the local approximation works if in- 
stead of a massive particle one has an extended 
object, small in comparison to its host galaxy? 

Textbook problems are designed in general to yield one 
correct answer, the above ideas for problems are rather 
vague but this is on purpose. The reason is twofold. On 
one hand, to promote in students a spirit of research by 
setting an approximate physical model and to look for the 
required data and "tools" to solve it; some of them can 
be found in the references. On the other hand, no single 
definite answer can be given. A feature proper of the 
way physics evolves toward describing and understanding 
nature. 



APPENDIX A: 



ASTRONOMICAL AND N-BODY 
UNITS 



Several quantities in astronomy are so large in compar- 
ison to common "terrestial" values, that special units are 
used. Table H| lists some of these and their equivalences 
in physical units. 

In the MKS system of units the Gravitational constant 
is G — 6.67 x 10 _11 m 3 kg -1 s~ 2 . A natural system 
of units for gravitational interactions is that where the 
gravitational constant is set to G = 1; in the same way 
as for quantum systems Planck's constant is usually set 
to h — 1. On dimensional grounds [G] — u^ui/u m ; where 
u m , ui, and u v correspond, respectively, to units of mass, 
length and velocity. 



G = 4.3007 x 10 



-3]ELP£ = 4.4984 xHT 3 - 

s 2 M m 



pc J 



MyrM G 



The transformation of G using length units such as kpc 
or Mpc (10 6 pc) is direct. Choosing ui and u m the unit 
of velocity and of time Ut, under an appropriate G value, 
are 



Glim 

ui 



Gu r< 



In this way the transformation from iV-body units, 
where G — M = R = 1, to physical ones can be made. 
Table UT) lists some values for different choices of ui and 
u m , and the resulting units of u v and ut- The entries 
correspond to using the approximate size and mass of a 
globular cluster, a disk of a spiral galaxy, and of a cluster 
of galaxies, respectively, as units ui and u m . 

In the standardized gravitational TV-body units&SS the 
total energy of a system is E = —1/4. This follows from 
the virial theorem (2T + W = 0), where 



2 R 



_ W _ GM 
~ T ~~ ~AR' 



Here R is strictly what is called the virial radius of the 
system; that does not necessarily coincides with the total 
extent of the stellar system, but is a very good approxi- 
mation. The potential energy of a Plummer model is 



p(r)ip(r)4:irr 2 dr = — 



3tt GM 2 
32 a 



Thus the total energy is E — -(3?rGM 2 )/(64a). In N- 
body units this leads to a value of the Plummer scale- 
radius of a — 37r/16. 
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